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O ' Abstract 

i I herence length (CL) effects and gluon shadowing (GS) we study electroproduction of 

charmonia off nuclei. In contrast to production of light vector mesons (p°, ( I )0 ) when at 
^ . small and medium energies CT and the onset of CL effects are not easily separated, in 

production of charmonia CT effects dominate. We found rather large CT effects in the 
range of Q 2 < 20 GeV 2 . They are stronger at low than at high energies and can be easily 
identified by the planned future experiments. Our parameter-free calculations explain 
well the NMC data for variation with photon energy of the Sn/C ratio of nuclear trans- 
parencies. We provide predictions for incoherent and coherent production of charmonia 
for future measurements. 
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Within a light-cone QCD formalism which incorporates color transparency (CT), co- 



1 



1 Introduction: space-time pattern of charmonium pro- 
duction 



The dynamics of charmonium production has been a hot topic evolved intensively during al- 
most three last decades. Discovery of J/^f in 1973 confirmed the idea of charm quark and 
gave a basis for its further investigations, which were affected also by new experiments sit- 
uated on more developed accelerators using more and more powerful electronics. Later, at 
the beginning of 90's the experiments with relativistic heavy-ion collisions stimulated the 
enhanced interest about charmonium suppression as a possible indication of the quark-gluon 
plasma formation. That fact has given a further motivation to continue in investigation of 
space-time pattern of charmonium production and opened new possibilities to analyze various 
corresponding phenomena. 

One of the fundamental phenomenon coming from QCD is color transparency (CT) studied 
intensively almost two last decades. This phenomenon can be treated either in the hadronic or 
in the quark basis. The former approach leads to Gribov's inelastic corrections ||, the latter 
one manifests itself as a result of color screening |3|, |3J]. Although these two approaches are 
complementary, the quark-gluon interpretation is more intuitive and straightforward. Colorless 
hadrons can interact only because color is distributed inside them. If the hadron transverse 
size r tends to zero then interaction cross section a(r) vanishes as r 2 ||. As a result the 
nuclear medium is more transparent for smaller transverse size of the hadron. Besides, this fact 
naturally explains the correlation between the cross sections of hadrons and their sizes || || 0j. 

Investigation of diffractive electroproduction of vector mesons off nuclei is very effective and 
sensitive for study of CT. A photon of high virtuality Q 2 is expected to produce a pair with 
a small ~ l/Q 2 transverse separation^. Then CT manifests itself as a vanishing absorption of 
the small sized colorless qq wave packet during propagation through the nucleus. Dynamical 
evolution of the small size qq pair to a normal size vector meson is controlled by the time scale 
called formation time. Due to uncertainty principle, one needs time interval to resolve different 
levels V (the ground state) or V (the next excited state) in the final state. In the rest frame 
of the nucleus this formation time is Lorentz dilated, 

*/ = — T ' 1 

TTly — Uly 

where v is the photon energy. A rigorous quantum-mechanical description of the pair evolution 
was suggested in || and is based on the light-cone Green function technique. A complementary 
description of the same process in the hadronic basis is presented in ||. 

Another phenomenon known to cause nuclear suppression is the effect of quantum coherence. 
It results from destructive interference of the amplitudes for which the interaction takes place on 
different bound nucleons. It reflects the distance from the absorption point when the pointlike 
photon becomes the hadronlike qq pair. This may be also interpreted as a lifetime of qq 

1 For production of light vector mesons (p°, <&°) very asymmetric pairs can be possible when the q or q carry 
almost the whole photon momentum. As a results the qq pair can have a large separation, see Sect. ^| and 
Eq. JIt]). Not so for production of charmonia where one can use the nonrelativistic approximation, a = 0.5, 
with rather high accuracy. 
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fluctuation providing the time scale which controls shadowing. Again, it can be estimated 
relying on the uncertainty principle and Lorentz time dilation as, 



tr 



2v 



Q 2 + m 2 v 



(2) 



It is usually called coherence time, but we also will use the term coherence length (CL), since 
light-cone kinematics is assumed, l c = t c (similarly, for formation length If = tf). CL is 
related to the longitudinal momentum transfer q c = l/l c in 7* iV — > V N, which controls the 
interference of the production amplitudes from different nucleons. 

Exclusive production of vector mesons at high energies is controlled by the sma\\-x Bj physics, 



and gluon shadowing becomes an important phenomenon [ ID| . It was shown in |Tl| that for 
electroproduction of charmonia off nuclei the gluon shadowing starts to be important at c.m.s. 
energy y/s > 30 — 60 GeV in dependence on nuclear target and Q 2 . Although the gluon 
shadowing within the kinematic range important for investigation of CT and discussed in the 
present paper is quite small we include it in all the calculations. 

In electroproduction of vector mesons off nuclei one needs to disentangle CT (absorption) 
and CL (shadowing) as the two sources of nuclear suppression. Detailed analysis of the CT 
and CL effects in electroproduction of vector mesons off nuclei showed that one can easily 
identify the difference of the nuclear suppression corresponding to absorption and shadowing 
in the two limiting cases for the example of vector dominance model (VDM). 
i.) The limit of Z c , shorter than the mean internucleon spacing ~ 2fm. In this case only 
final state absorption matters. The ratio of the quasielastic (or incoherent) 7* A — > V X and 
7* iV — > V X cross sections, usually called nuclear transparency, reads H, 
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Here z is the longitudinal coordinate and b is the impact parameter of the production point of 
vector meson. In (§) pA^b, z) is the nuclear density and a^ N is the inelastic VN cross section, 
ii.) In the limit of long l c expression for nuclear transparency takes a different form, 



Tr mc 



d 2 bT A (b) exp -a^T A (b) 



(4) 



where we assume a^ N aY n N for the sake of simplicity. T A {b) is the nuclear thickness function 

00 

T A (b) = J dzp A (b,z) . (5) 

—00 

The exact expression which interpolates between the two regimes ([3]) and (f|) can be found in 



0- 



The problem of separation of CT and CL effects arises especially in production of light vector 
mesons (p°, $°) [10|. In this case the coherence length is larger or comparable with the formation 
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one starting from the photoproduction limit up to Q 2 ~ l-i-2 GeV 2 . In charmonium production, 
however, there is a strong inequality If > l c independently of Q 2 and v. It leads to a different 
scenario of CT-CL mixing as compared with production of light vector mesons. That fact gives a 
motivation for separated study of J/\l/ production presented in this paper using light-cone dipole 



approach generalized for the case of a finite coherence length and developed in [fTOj] . Another 
reason is supported by the recent paper where charmonium production was calculated in 
the approximation of long coherence length l c 3> Ra using realistic charmonia wave functions 
and corrections for finite values of l c . It gives very interesting possibility to compare the 
predictions of the present paper with the results obtained from Jll]| for enhancement of reliability 
of theoretical predictions as a realistic basis for planned future experiments of electron-nucleus 
collisions. 

The paper is organized as follows. In Sect. ^| we present a short review of the light-cone 
(LC) approach to diffractive electroproduction of vector mesons in the rest frame of the nucleon 
target. Here we also present the individual ingredients contained in the production amplitude: 
(i) the dipole cross section characterizing the universal interaction cross section for a color- 
less quark-antiquark dipole and a nucleon. (ii) The LC wave function for a quark-antiquark 
fluctuation of the virtual photon, (iii) The LC wave function of charmonia. 

As the first test of the model we calculate in Sect. | the cross section of elastic electropro- 
duction of J/\l/ off a nucleon target. These parameter-free calculations reproduce both energy 
and Q 2 dependence remarkably well, including the absolute normalization. 

Sect. [| is devoted to incoherent production of J/\I / off nuclei. Numerical calculations are 
compared with the available NMC data for variation with photon energy of the Sn/C ratio 
of nuclear transparencies. We find a different scenario of an interplay between coherence and 
formation length effects from one occurred in production of light vector mesons. Because a 
variation of l c with Q 2 can mimic CT at medium and low energies, one can map experimental 
events in Q 2 and v in such a way as to keep l c = const. The LC dipole formalism predicts 
rather large effect of CT in the range of Q 2 < 20 GeV 2 . This fact makes it feasible to find a 
clear signal of CT effects also in exclusive production of J/\l/ in the planned future experiments. 

Coherent production of vector mesons off nuclei leaving the nucleus intact is studied in 
Sect. |^. The detailed calculations show that the effect of CT on the Q 2 dependence of nuclear 
transparency at l c = const is weaker than in the case of incoherent production and is difficult 
to be detected at low energies since the cross section is small. 

Although the gluon shadowing starts to manifest itself at \fs > 30 — 60 GeV and is not very 
significant in the energy range important for searching for CT effects in the present paper we 
include it in all calculations. 

The results of the paper are summarized and discussed in Sect. [| An optimistic prognosis 
for discovery of CT in electroproduction of charmonia is made for the future experiments. 
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2 A short review of the light-cone dipole phenomenol- 
ogy for elastic electroproduction of charmonia 7*AT — > 
J/*ff N 

The light-cone (LC) dipole approach for elastic electroproduction 7*iV — > V N was already 
used in papers for study of exclusive photo- and electroproduction of charmonia and in [[10] 
for elastic virtual photoproduction of light vector mesons p° and $°. Therefore, we present only 
a short review of this LC phenomenology with the main emphasis to elastic electroproduction 
of charmonia. Here a diffractive process is treated as elastic scattering of a qq fluctuation 
(cc fluctuation for the case of charmonium production) of the incident particle. The elastic 
amplitude is given by convolution of the universal flavor independent dipole cross section for 
the qq interaction with a nucleon, a qq , || and the initial and final wave functions. For the 
exclusive photo- or electroproduction of charmonia 7*iV — > J/ty N the forward production 
amplitude is represented in the quantum-mechanical form 

i 

M^n^j/v N {s, Q 2 ) = (J/V\a£(r, s)| 7 *) = J da J d 2 r *} /v (r, a) a qq (f, s) # 7 *(f, a, Q 2 ) (6) 
with the normalization 
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In order to calculate the photoproduction amplitude one needs to know the following in- 
gredients of Eq. (^): (i) the dipole cross section a qq (r,s) which depends on the qq transverse 
separation f and the cm. energy squared s. (ii) The light-cone (LC) wave function of the 
photon \l/ 7 * (r*, a, Q 2 ) which also depends on the photon virtuality Q 2 and the relative share a 
of the photon momentum carried by the quark, (iii) The LC wave function j/^(r, a) of J . 

Note that in the LC formalism the photon and meson wave functions contain also higher 
Fock states \qq), \qqG), \qq2G), etc. The effects of higher Fock states are implicitly incorporated 
into the energy dependence of the dipole cross section <y qq (r, s) as is given in Eq. (||). 

The dipole cross section a qq (r,s) represents the interaction of a qq dipole of transverse 
separation f with a nucleon ||. It is a flavor independent universal function of r and energy 
and allows to describe in a uniform way various high energy processes. It is known to vanish 
quadratically a qq (r,s) oc r 2 as r — > due to color screening (CT property) and cannot be 
predicted reliably because of poorly known higher order pQCD corrections and nonperturbative 
effects. Detailed discussion about the dipole cross section u qq (r, s) with emphasis to production 
of light vector mesons is presented in Jl0| . In electroproduction of charmonia the corresponding 



transverse separations of cc-dipole reach the values < 0.4 fm (semiperturbative region). It means 
that nonperturbative effects are sufficiently smaller as compared with light vector mesons. 
Similarly, the relativistic correstions are also small enough to use safely the nonrelativistic limit 
a = 0.5 with rather high accuracy 0. 

There are two popular parameterizations of cr qq {r, s). The first one suggested in [TJ[] reflects 



the fact that at small separations the dipole cross section should be a function of r and XBj ~ 
l/(r 2 s) to reproduce Bjorken scaling. It well describes data for DIS at small x and medium 
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and high Q 2 . However, at small Q 2 it cannot be correct since it predicts energy independent 
hadronic cross sections. Besides, Xsj is not any more a proper variable at small Q 2 and should 



be replaced by energy. This defect is removed by the second parameterization suggested in [15 



which is similar to one in [jT3] , but contains an explicit dependence on energy. It is valid down 
to the limit of real photoproduction. Since we want to study CT effects starting from Q 2 = 0, 
we choose the second parametrization, which has the following form : 



cr, 



qq 



(r,s) = a (s) 
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,Sq, 



(8) 
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Here (r 2 h ) = 0.44 fm 2 is the mean pion charge radius squared; Sq = 1000 GeV 2 . The cross 



section cr^(s) was fitted to data in [16, 17], 



aZ{s) = 23.6 



0.079 



mb 



(11) 



The dipole cross section Eqs.(|8]) - ( |i~T| ) provides the imaginary part of the elastic amplitude. 
It is known, however, that the energy dependence of the total cross section generates also a 



real part [13], 



1 — % — 



d 



(12) 



is rather steep at small r leading to 



2 ain(s) 

The energy dependence of the dipole cross section Eq. 
a large real part which should not be neglected. For instance, the photoproduction amplitude 
of 7iV — > J/tyN rises oc s a2 and the real-to- imaginary part ratio is over 30%. 

Although the form of Eq. @ successfully describes data for DIS at small x only up to 
Q 2 fa 10 GeV 2 we prefer this parameterization for study of charmonium electroproduction. 
The reason is that we want to study CT effects predominantly in the range Q 2 < 20 GeV 2 
and in addition parameterization Eq. (^j) describes the transition toward photoproduction limit 
better than parameterization presented in ||14|| . Besides, in the paper El was shown study- 



ing electroproduction of charmonia off nucleons that the difference between predictions using 
both parameterizations [13] and Eq. ([8]) is rather small and can be taken as a measure of the 
theoretical uncertainty. 



The perturbative distribution amplitude ("wave function") of the qq (cc for J/^/ production) 
Fock component of the photon has the following form for transversely (T) and longitudinally 
(L) polarized photons p|, ^ ^TJ, 



*£ L (r» = Z q xO T > L X K (er) (13) 
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where \ and x are the spinors of the quark and antiquark, respectively; Z q is the quark charge, 
Z q = Z c = 2/3 for J/^f production; Nq = 3 is the number of colors. K (er) is a modified Bessel 
function with 

e 2 = a(l-a)Q 2 + m 2 c , (14) 

where m c = 1.5 GeV is mass of the c quark, and a is the fraction of the LC momentum of the 
photon carried by the quark. The operators T)L read, 



O 



r 



m c a ■ e + i (1 — 2a) (<r • n) (e • V r ) + (a x e) ■ V r 



L = 2Qa{\ - a) 



[a ■ n) 



(15) 
(16) 



Here V r acts on transverse coordinate r; e is the polarization vector of the photon and n is a 
unit vector parallel to the photon momentum. 

In general, the transverse qq separation is controlled by the distribution amplitude Eq. ( [T3|) 
with the mean value, 

(r) ~ - = 1 = . (17) 

For production of light vector meson very asymmetric qq pairs with a or (1 — a) "%jVP? q jQ 2 
become possible. Consequently, the mean transverse separation (r) ~ l/m q becomes huge since 
one must use current quark masses within pQCD. However, that is not the case in charmonium 
production because of a large quark mass m c = 1.5 GeV. Therefore, we are out of the problem 
how to include nonperturbative interaction effects between c and c because they are rather 
small. Despite of this fact for completeness we include these nonperturbative interaction effects 
in all calculations to avoid although small but supplementary uncertainties in predictions. We 
take from |T5] the corresponding phenomenology including the interaction between c and c 



based on the light-cone Green function approach. 

The Green function G qq (zi, t\\ z 2 , r 2 ) describes the propagation of an interacting qq pair (cc 
pair for the case of J/^> production) between points with longitudinal coordinates z\ and z 2 and 
with initial and final separations f\ and r*2. This Green function satisfies the two-dimensional 
Schrodinger equation, 



d 

dzo 



Gqqiz^n; z 2 ,f 2 ) 



e 2 - A T 



2 va, (1 — a] 



+ V qq (z 2 ,f 2 ,a) 



G qq (zi,fi; z 2 ,f 2 ) 



Here v is the photon energy. The Laplacian A r acts on the coordinate r. 

The imaginary part of the LC potential V qq (z 2 , f 2 , a) in ([Tj|) is responsible for attenuation 
of the qq in the medium, while the real part represents the interaction between the q and q. 
This potential is supposed to provide the correct LC wave functions of vector mesons. For the 
sake of simplicity we use the oscillator form of the potential, 



Re VgJz 2 , r 2 , a) 



a 4 (a) r 2 2 
2 v a(l — a) ' 



(19) 



which leads to a Gaussian r-dependence of the LC wave function of the meson ground state. 
The shape of the function a(a) will be discussed below. 
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In this case equation (|i"8"D has an analytical solution, the harmonic oscillator Green function 



22 
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x exp 
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where Az = z 2 — z\ and 



a 2 (a) 



(20) 



(21) 



v a(l — a) 

The boundary condition is Gq q (zi,fi; z 2 ,r 2 )\ z . 2=Zl = b 2 (fi — r 2 ). 

The probability amplitude to find the qq fluctuation of a photon at the point z 2 with 
separation f is given by an integral over the point Z\ where the qq is created by the photon 
with initial separation zero, 



*£ L (r,a)- /Z " V "' 



47r v a(l — a) 



dzi (x T,L x) Gq q (zi, n;z 2 , r) 



ri=0 



(22) 



The operators T)L are defined in Eqs. ([15]) and (0). Here they act on the coordinate T\. 
If we write the transverse part as 



X0 T X = A + B- V r 



then the distribution functions read, 



a) = Zqy/^ A%(e,r,\) + 5$ 1 (e,r,A) 



VL(r, a) = 2 ZqiJa^Q a(l - a) x a • n X $o(e, r, A) 



where 
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2a 2 (a) 



The functions $ ,i i n Eqs. ( p4| ) and (pq) are defined as 
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(23) 

(24) 
(25) 

(26) 

(27) 
(28) 



Note that the q — q interaction enters Eqs. and ( Pq) via the parameter A defined in 
In the limit of vanishing interaction A — * (i.e. Q 2 — > oo, a is fixed, a ^ or 1) Eqs. fl24| ) - (p5|) 
produce the perturbative expressions of Eq. (O)- As was mentioned above, for charmonium 
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production nonperturbative interaction effects are quite weak. Consequently, the parameter A 
is rather small due to a large quark mass m c — 1.5 GeV : 

Q z + 4 

With the choice a 2 (a) oc a(l — a) the end-point behavior of the mean square interquark 
separation (r 2 ) oc l/a(l — a) contradicts the idea of confinement. Following [|15j we fix this 
problem via a simple modification of the LC potential, 

a 2 (a) = Oq + 4a 2 a(l — a) . (30) 



The parameters ao and ai were adjusted in |15| to data on total photoabsorption cross section 



23, |2l|, diffractive photon dissociation and shadowing in nuclear photoabsorption reaction. 
The results of our calculations vary within only 1% when ao and a\ satisfy the relation, 

a 2 = V 115 (0.112) 2 GeV 2 

a\ = (1 — l>) L15 (0.165) 2 GeV 2 , (31) 

where v takes any value < v < 1. In view of this insensitivity of the observables we fix the 
parameters at v = 1/2. We checked that this choice does not affect our results beyond a few 
percent uncertainty. 



The last ingredient in elastic production amplitude @ is charmonium wave function. We 
use a popular prescription |25| applying the Lorentz boost to the rest frame wave function 
assumed to be Gaussian which leads to radial parts of transversely and longitudinally polarized 
mesons in the form, 

a(l — a) r 2 

^{r,a) = U~>~a[L -a)f(a)exp 



2B 2 



(32) 



with a normalization defined below, and 

/(«) = 



exp 



m 2 R 2 



(33) 



2a(l - a) 

with the parameters from [25||, R = 0.183 fm and m c = 1.3 GeV. A detailed analysis of various 
problems in this relativization procedure |27j leads to the same form d32|). 

We assume that the distribution amplitude of cc fluctuations for J/\& and for the photon 
have a similar structure ||26|| . Then in analogy to Eqs. 



[r,a) 



(A + B-V)<S> T J/q ,(r } a) 



(34) 

2m J/qi a(l-a)(x^-nx)^j/^(r,a) . (35) 
Correspondingly, the normalization conditions for the transverse and longitudinal charmo- 



nium wave functions read, 
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(36) 
(37) 
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3 Electroproduction of J/ty on a nucleon, comparison 
with data 



In this section we verify first the LC approach by comparing with data for nucleon target. 
The forward production amplitude 7* iV — > J/ty N for transverse and longitudinal photons and 
charmonium is calculated using the nonperturbative photon Eqs. (p4|), (|25|) and vector meson 
Eqs. (|34]), ( |35| ) wave functions, and has the following form, 



1 

M^ N ^j/^ N (s,Q 2 )\ t=o = N c Z cy fa^ J d 2 ra qq (f,s) J da{m 2 c <S> (e,r, \)^ 



+ \a 2 + (l-a) 2 



$i(e,F,A)-V r #J/*(r,a)} 



(3f 



s,Q 2 ) t=Q = 4:N C Z cy /a^mj/yQ Jd 2 ra qq (r,s) 



x / da a (1 - at) <E> (e, f, X)^j /9 (f, a) 


T,L 



(39) 



These amplitudes are normalized as | Ai T ' L \ 2 = IQtc da^j jdt . The real part of the amplitude 



t=o 



is included according to the prescription described in the previous section. We calculate the 
cross sections a = a T + ea L assuming that the photon polarization is e = 1. 

Now we can check the absolute value of the production cross section by comparing with 
data for elastic charmonium electroproduction 7*^9 — > J/^fr p. Unfortunately, data are available 
only for the cross section integrated over t, 



a 



( 7 *iV -> VN) 



\M 



T,L\2 



7* AT 



(40) 



J/ty p. We use the experimental value [p8 



where B is the slope parameter in reaction 7*p 
B = 4.7 GeV" 2 . 

Our predictions are plotted in Fig. |l] together with the data on the Q 2 + m 2 ,^ dependence 



of the cross section from HI and ZEUS JSIJ. 

As the second test of our approach is a description of the low Q 2 region as well. As we 
discussed in the previous section despite of smallness of effects of nonperturbative interaction 
between the c and c, we include them into calculations. Comparison of the model with data 
PT| , |31| |3"3H for the energy dependence of the cross section of real J /\& photoproduction is 
presented in Fig. ^. 

The normalization of the cross section and its energy and Q 2 dependence are remarkably well 
reproduced in Figs. [1] - 12|. This is an important achievement since the absolute normalization is 
usually much more difficult to reproduce the production cross sections than nuclear effects. For 
instance, the similar, but simplified calculations in || underestimate the J/\I/ photoproduction 
cross section on protons by an order of magnitude. 

As a cross-check for the choice of the J/\l/ wave function in Eqs. (|3"2D and (|3"3"|) we also 



calculated the total J/\l/-nucleon cross section, which was already estimated in |13[ using the 
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Figure 1: Q 2 + m 2 ^- dependence of the integrated cross section for the 
reactions ^* p — > J/ty p. The model calculations are compared with HI 
m and ZEUS data at energy W = 90 GeV. 



charmonium wave functions calculated with several realistic qq potentials. The J/\l/-nucleon 
total cross section has the form, 



lot 



N, 



c 



d r da <m 



+ [a 2 + (1 - a) 



(41) 



We calculated <Jtot S ' N with the charmonium wave function in the form (^) with the parameters 
described in the previous section. For the dipole cross section we adopt the KST parame- 
terization (H) which is designed to describe low-Q 2 data. Then, at yfs = 10 GeV we obtain 
4.2 mb which is in a good agreement with <J^ N 



J/^N 
&tot 



3.6 mb evaluated in [13] 



4 Incoherent production of charmonia off nuclei 
4.1 the LC Green function formalism 

In this section we present a short review of the LC Green function formalism for incoherent 
production of an arbitrary vector meson. For the case of charmonium production one should 
replace V — > J/\I/ and qq — > cc. In general, the diffractive incoherent (quasielastic) production 
of vector mesons off nuclei, 7* A — > VX, is associated with a break up of the nucleus, but 
without production of new particles. In another words, one sums over all final states of the 
target nucleus except those which contain particle (pion) creation. The observable usually 
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Figure 2: Energy dependence of the real photoproduction cross section 
on a nucleon, 7 p — > J/ty p. Our results are compared with data from the 
fixed target E401 j^Tjj, E516 g2j, and collider HERA HI 
experiments. 



and ZEUS 



studied experimentally is nuclear transparency denned as 
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(42) 



The t-slope of the differential quasielastic cross section is the same as on a nucleon target. 
Therefore, instead of integrated cross sections one can also use nuclear transparency expressed 
via the forward differential cross sections Eq. (^), 
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j*A^VX 



[s,Q 2 



M 1 *n->vn(s } Q' 2 



(43) 



In the LC Green function approach [|T^] the physical photon |y) is decomposed into different 
Fock states, namely, the bare photon |7*)o, | <?<?), \qqG), etc. As we mentioned above the higher 
Fock states containing gluons describe the energy dependence of the photoproduction reaction 
on a nucleon. Besides, those Fock components also lead to gluon shadowing as far as nuclear 
effects are concerned. However, these fluctuations are heavier and have a shorter coherence 
time (lifetime) than the lowest \qq) state. Therefore, at medium energies only \qq) fluctuations 
of the photon matter. Consequently, gluon shadowing related to the higher Fock states will 
be dominated at high energies. Detailed description and calculation of gluon shadowing for 
the case of vector meson production off nuclei is presented in [110, O. Although the gluon 
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shadowing is quite small in the kinematic range for study of CT effects suggested in the present 
paper we include it in all the calculations. 

Propagation of an interacting qq pair in a nuclear medium is also described by the Green 
function satisfying the evolution Eq. fll8|) . However, the potential in this case acquires an 
imaginary part which represents absorption in the medium (see (||) for notations), 



ImVq q (z 2 ,r,a) 



PaQ>, z 2 ) 



(44) 



The evolution equation ( |T8D with the potential Vq q (z 2 , f 2 , a) containing this imaginary part 
was used in f34] , |35| , and nuclear shadowing in deep-inelastic scattering was calculated in good 
agreement with data. 

The analytical solution of Eq. ( |20| ) is only known for the harmonic oscillator potential V(r) oc 
r 2 . To keep the calculations reasonably simple we are forced to use the dipole approximation 



a qq( r ' S ) = C(s)r 2 



(45) 



which allows to obtain the Green function in an analytical form. 

The energy dependent factor C(s) is adjusted by demanding that calculations employing 
the approximation Eq. flSSf ) reproduce correctly the results based on the realistic cross section 
in the limit l c 3> R A (the so called "frozen" approximation) when the Green function takes the 
simple form : 



Gq q {zi, ri; z 2 , f 2 ) =^6{i\- f 2 ) exp 



-(Jqqin) I dzp A (b,z) 



(46) 



where the dependence of the Green function on impact parameter is dropped. Thus, for inco- 
herent production of vector mesons the factor C(s) is fixed by the relation, 



Jd 2 bT A (b) 



J d 2 rr 2 exp 



-\C T ^s)r 2 T A {b) 



J da^* v (r,a) ^i L (r,a 



J d 2 r r 2 J da^y(f,a) ^ 



T.Li 



[r,a, 



Jd 2 bT A (b) 



J d 2 r <Jq q {r, s) exp 



\aqq{r,s)T A {b) 



J d 2 r aq q (r, s) J da ^(r, a) ^ L (f, a 



(47) 



To take advantage of the analytical form of the Green function which is known only for the 
LC potential Eq. ([0]) with a constant nuclear density, we use the approximation p A (b,z) = 
Po Q{R A — b 2 — z 2 ). Therefore we have to use this form for Eq. ([|?]) as well. The value of the 
mean nuclear density po has been determined using the relation, 



d 2 b 



1 - expi -co p jR 2 A - b 2 



d 2 b 



1 — exp 



00 



T(b) 



(4? 



where the nuclear thickness function T A (b) is calculated with the realistic Wood-Saxon form of 
the nuclear density. The value of p turns out to be practically independent of the cross section 
o"o in the range from 1 to 50 mb. 
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With the potential Eqs. (44) - (|45| ) the solution of Eq. ([Tj|) has the same form as Eq. (|20|) , 
except that one should replace u =>■ Q, where 



n 



J a 4 (a) — i paQ>, z) v a (1 — a) C(s) 



v ail 



a) 



(49) 



As we discussed in [[HJ the value of l c can distinguish different regimes of vector meson 
production. 

(i) The CL is much shorter than the mean nucleon spacing in a nucleus (l c — > 0). In this case 
G(z 2 ,f 2 ; Zi,ri) — > 5(z 2 — Z\). Correspondingly, the formation time of the meson wave function 
is very short as well as given in Eq. ([!]). For light vector mesons If ~ l c and since formation 
and coherence lengths are proportional to photon energy both must be short. Consequently, 
nuclear transparency is given by the simple formula Eq. corresponding to the Glauber 
approximation. 

(ii) In production of charmonia and other heavy flavor quarkonii the intermediate case 
l c — ► 0, but If ~ Ra can be realized. Then the formation of the meson wave function is 
described by the Green function and the numerator of the nuclear transparency ratio Eq. (ff~3|) 
has the form M, 



M 



-y*A^VX 



s,Q 2 



dzp A (b,z) Ftfaz) 



where 



F 1 (b,z) 



da / d 2 n d 2 r 2 ^* v {f 2} a) G(z\ r 2 ; z, n) a qq {r 1} s) * 7 * (n, a) 



(50) 



(51) 



(iii) In the high energy limit l c 3> Ra (in fact, it is more correct to compare with the mean 
free path of the qq in a nuclear medium if the latter is shorter than the nuclear radius). In 
this case G(z 2 ,r 2 ; 2a,r0 — > 5{r 2 — fx), i.e. all fluctuations of the transverse qq separation are 
"frozen" by Lorentz time dilation. Then, the numerator on the r.h.s. of Eq. (H3) takes the form 



M 



■y*A^VX 



d 2 bT A (b) 



d 2 r / da 
Jo 



(52) 



x tyy(r,a)<7q q (r,s) exp 



--(Tq q (r, s)T A (b) 



* r (f, a, Q 2 



In this case the qq attenuates with a constant absorption cross section like in the Glauber 
model, except that the whole exponential is averaged rather than just the cross section in 
the exponent. The difference between the results of the two prescriptions are the well known 
inelastic corrections of Gribov ||. 

(iv) This regime reflects the general case when there is no restrictions for either l c or If. 



The corresponding theoretical tool has been developed for the first time only recently in [10 



and applied to electroproduction of light vector mesons at medium and high energies. Even 
within the VDM the Glauber model expression interpolating between the limiting cases of low 
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[(i), (ii)] and high [(iii)] energies has been derived only recently [12] as well. In this general case 
the incoherent photoproduction amplitude is represented as a sum of two terms [|36|. 



M 



~/*A->VX\ 



d 2 b / dzp A (b,z) F x (b, z) - F 2 (b, z 



(53) 



The first term F\ (6, z) introduced above in Eq. ( |51"D alone would correspond to the short l c limit 
(ii). The second term F 2 (b, z) in fl53|) corresponds to the situation when the incident photon 
produces a qq pair diffractively and coherently at the point Z\ prior to incoherent quasielastic 
scattering at point z. The LC Green functions describe the evolution of the qq over the distance 
from Z\ to z and further on, up to the formation of the meson wave function. Correspondingly, 
this term has the form, 



F 2 (b,z) 



_L 

dz\ pA(b, z\) J "da j d 2 rid 2 r 2 d 2 r^ v (r 2 ,a) 



x G{z -> oo, r 2 ; z, r) (Tg q (f, s) G(z, f; z 1 ,r 1 ) <7 ft (ri, s) # 7 * (r x , a) . 



(54) 



Eq. ( |5*B"D correctly reproduces the limits (i) - (iii). At l c — > the second term F 2 (b, z) 
vanishes because of strong oscillations, and Eq. (|53|) reproduces the Glauber expression Eq. (|3]). 
At l c 3> Ra the phase shift in the Green functions can be neglected and they acquire the simple 
form G(z 2 ,f 2 ;zi,fi) — > 5(f 2 — r*i). In this case the integration over longitudinal coordinates 
in Eqs. (fy]) and ( |54"1) can be performed explicitly and the asymptotic expression Eq. (|52"D is 
recovered as well. 



4.2 Comparison with data for incoherent production of J 



Exclusive incoherent electroproduction of vector mesons off nuclei has been suggested in |J7| to 
be very convenient for investigation of CT. Increasing the photon virtuality Q 2 one squeezes the 
produced qq wave packet. Such a small colorless system propagates through the nucleus with 
little attenuation, provided that the energy is sufficiently high (// 3> Ra) so the fluctuations of 
the qq separation are frozen during propagation. Consequently, a rise of nuclear transparency 
Tr l ^ c (Q 2 ) with Q 2 should give a signal for CT. Indeed, such a rise was observed in the E665 
experiment |38| at Fermilab for exclusive production of p° mesons off nuclei what has been 



claimed as manifestation of CT. 

However, the effect of coherence length |$9|, leads also to a rise of Tr l ^ c {Q 2 ) with Q 2 



and so can imitate CT effects. This happens when the coherence length varies from long to 
short (see Eq. (|2])) compared to the nuclear size and the length of the path in nuclear matter 
becomes shorter. Consequently, the vector meson (or qq) attenuates less in nuclear medium. 
This happens when Q 2 increases at fixed v. Therefore one should carefully disentangle these 
two phenomena. 

Unfortunately the data on electroproduction of charmonia are very scanty so far. There are 
only data from the NMC experiment 0] concerning energy dependence of the ratio of nuclear 



transparencies Tr 1 ^ and Tr^ c for incoherent production of J/\l/ at Q 2 = 0. The corresponding 
photon energy varies from 60 to 210 GeV. It allows to study the transition from medium long 
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Figure 3: Energy dependence of the ratio of nuclear transparencies Tr 



Sn 



and Tr<7 vs. experimental points taken from the NMC experiment [40 



Solid and dashed curves show our results using the LC Green function 
approach in general case with no restriction for either l c or If Eq. ( |53"D 
and in the limit of l c — > Eq. (150), respectively. 



to long coherence length, which varies from 2.4 to 8.5 fm. For long l c > 8.5 fm the "frozen" 
approximation can be used with high accuracy. In this case nuclear transparency Tr™ c of 
incoherent (quasielastic) J/\l/ production can be calculated using Eq. (|53D and the simplified 
"frozen" approximation Eqs. (|4TJ) - For medium long coherence length one can not use the 
"frozen" approximation and fluctuations of the size of the qq pair become important. Because of 
a strong inequality l c < If for charmonium production CT effects are expected to be dominant 
at small and moderate energies. Consequently, they should lead to a rise with energy of Tr™ c . 
Such a scenario is depicted in Fig. ||| by solid and dashed curves. Dashed curve show our results 
using the LC Green function approach in the limit of short coherence length l c — > Eq. (|50|). 
The solid one includes in addition also CL effects. Thus, the effect of coherence length manifest 
itself as a separation between the solid and dashed curves. A rise with energy of the ratio 
Tr™^ /Tr™ c at small and medium energy is a net manifestation of CT. It follows from a rise of 
formation time, see Eq. (JI|). At larger energies when CL effects become also important the ratio 
Tr^/Tr™ starts to be gradually smaller^. Unfortunately, the NMC data have quite large error 
bars and therefore give only an indication for such a behavior. More accurate data are needed 
for exploratory study of CT and CL effects. Charmonium real photoproduction off nuclei at 

2 In energy dependence of nuclear transparency at fixed Q 2 effect of the coherence follows from variation of 
the coherence length from small to large values compared to the nuclear size, see Eq. (||) 
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small and large energies is very sensitive for investigation of CT and CL effect separately. Not 
so for real photoproduction of light vector mesons when coherence and formation lengths are 
comparable and mixing of CT and CL effects exists already at small energies. 



Problem of separation of CT and CL effects was discussed in details in [IIJ with the main 
emphasis to production of light vector mesons where l c > If at Q 2 < 1 4- 2 GeV 2 . In this 
paper we present the results for charmonium production where a strong inequality l c < If in 
all discussed kinematic region leads to a different scenario of CT-CL mixing as compared with 
light vector mesons. Consequently, at fixed Q 2 and at small and medium energies the problem 
of CT-CL separation is not so arisen. In order to eliminate the effect of CL from the data on 
the Q 2 dependence of nuclear transparency one should bin the data in a such way which keeps 
l c = const |J. It means that one should vary simultaneously v and Q 2 maintaining the CL 
Eq. (H) constant, 

v = \lc {Q 2 + m 2 ^) . (55) 

In this case the Glauber model predicts a Q 2 independent nuclear transparency, and any rise 
with Q 2 would signal CT §. 

The LC Green function technique incorporates both the effects of coherence and formation. 
We performed calculations of Tr l ^ c (Q 2 ) at fixed l c starting from different minimal values of v, 
which correspond to real photoproduction in Eq. (|55|), 

Vmin = ~2^ c m J/* ' (56) 

The results for incoherent production of J/ty at v m in = 24.3, 121.7 and 487 GeV (Z c = 1,5 and 
20 fm) are presented in Fig. ^ for beryllium, iron and lead. We use the nonperturbative LC 
wave function of the photon with the parameters of the LC potential ao,i fixed in accordance 
with Eq. (|3T|) at v — 1/2. We use quark mass m c = 1.5 GeV. 

Although the predicted variation of nuclear transparency with Q 2 at fixed l c is less than for 



production of light vector mesons [10] it is still rather significant to be investigated experimen- 
tally even in the range of Q 2 < 20 GeV 2 . CT effects (the rise with Q 2 of nuclear transparency) 
are more pronounced at low than at high energies and can be easily identified by the planned 
future experiments. 

We also calculated the energy dependence of nuclear transparency at fixed Q 2 . The results 
for beryllium, iron and lead are shown in Fig. [| for different values of Q 2 . The interesting 
feature is the presence of a maximum of transparency at some energy which is much more 
evident than in production of light vector mesons flQfl . At small and moderate energies a 
strong rise of Tr™ c with energy is a manifestation of net CT effects as a result of a strong 
inequality l c < If. The existence of maxima of Tr™ results from the interplay of coherence and 
formation effects. Indeed, the formation length (FL) rises with energy leading to an increasing 
nuclear transparency. At some energy, however, the effect of CL which is shorter than the 
FL, is switched on leading to a growth of the path length of the qq in the nucleus, i.e. to a 
suppression of nuclear transparency. This also explains the unusual ordering of curves calculated 
for different values of l c as is depicted in Fig. £|. 
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5 Coherent production of J/ty 



First of all we present a short introduction to coherent production of vector mesons. One should 
replace V — > J/\& and qq — > cc when coherent production of charmonia is treated. In general, 
in coherent (elastic) electroproduction of a vector mesons the target nucleus remains intact, 
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Figure 4: Q 2 dependence of the nuclear transparency Tr^ lc for exclusive 
electroproduction of J/\& on nuclear targets 9 -Be, 56 Fe and 207 Pb (from 
top to bottom). The CL is fixed at l c — 1, 5 and 20 fm. 



18 



so all the vector mesons produced at different longitudinal coordinates and impact parameters 
add up coherently. This condition considerably simplifies the expressions for the production 
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Figure 5: Nuclear transparency for incoherent electroproduction 7* 'A — > 
J/* X as function of energy at Q 2 = 0, 5, 20 and 100 GeV 2 for beryllium, 
iron and lead. The solid curves and dashed curves for lead correspond to 
calculations with and without gluon shadowing, respectively. 
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cross sections. The integrated cross section has the form, 



-.coh 



a. 



coh 

-y*A^VA 



d 2 q 



d 2 be ig.b M coh A _^ vA ^ 



d 2 b\M^ VA {b) 



where 



M^ VA (b) 



dz p A (b,z) Fx(b,z) 



(57) 



(5? 



with the function F\(b, z) defined in (|51|). 

One should not use Eq. (f43D for nuclear transparency any more since the t-slopes of the 
differential cross sections for nucleon and nuclear targets are different and do not cancel in 
the ratio. Therefore, the nuclear transparency also includes the slope parameter B^*n for the 
process 7* N — > V N, 



Tr 



coh 



a 



coh 



16ir B 7 *n &a 



coh 



Aa N A\M rN ^ V N(s,Q 2 )\2 (59) 

The energy dependent factor C(s) in dipole cross section approximation Eq. fl45l) is adjusted 
in the limit l c 3> Ra in analogical way as for incoherent production of charmonia described 
in the previous section. However, in contrast to Eq. (]47|) the factor C(s) is fixed now by the 
following relation, 



fd 2 b 



J d 2 r J da V^(r, a) #5 L (r, a) \ l - exp 



-\C T r{s)r 2 T A {b) 



J d?r J da V* v (r, a) C t ,l(s) r 2 ^l L (f, a 



Jd 2 b 



J d 2 r J da *^(f, a) ^ L (f, a) \ l - exp 



-l°q q (r,s)T A (b) 



J d?r J da y* v (r, a) a qq (r, s) ^ L (f, a) 



(60) 



5.1 Predictions for coherent production of J/ty 

Unfortunately, there are no data yet on coherent electroproduction of charmonia therefore we 
present only predictions which can be later verify and tested in the future planned experiments. 

One can eliminate the effects of CL and single out the net CT effect in a way similar to 
what was suggested for incoherent reactions by selecting experimental events with l c = const. 
We calculated nuclear transparency for the coherent reaction 7* A —* J/^A at fixed values of 
l c . The results for l c = 1, 5 and 20 fm are depicted in Fig. |6| for several nuclei. We performed 
calculations of Tr c A oh with the slope B = 4.7 GeV -2 . The effect of a rise of Tr A h is not sufficiently 
large to be observable in the range Q 2 < 20GeV 2 . A wider range Q 2 < 100 GeV 2 and heavy 
nuclei gives higher chances for experimental investigation of CT. However, on the other hand, 
it encounters the problem of low yields at high Q 2 . 

Note that in contrast to incoherent production where nuclear transparency is expected to 
saturate as Tr l A c (Q 2 ) — ► 1 at large Q 2 , for the coherent process nuclear transparency reaches a 
higher limit, Tr A h (Q 2 ) -> A 1/3 . 

We also calculated nuclear transparency as function of energy at fixed Q 2 . The results for 
J/ty produced coherently off beryllium, iron and lead are depicted in Fig. [7| at Q 2 = 0, 5, 20 
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and 100 GeV 2 . Tr c A oh is very small at low energy, what of course does not mean that nuclear 
matter is not transparent, but the nuclear coherent cross section is suppressed by the nuclear 
form factor. Indeed, the longitudinal momentum transfer which is equal to the inverse CL, is 
large when the CL is short. However, at high energy l c 3> Ra and nuclear transparency nearly 
saturates (it decreases with v only due to the rising dipole cross section). The saturation level 
is higher at larger Q 2 which is a manifestation of CT. 
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Note that in all the calculations the effects of gluon shadowing are included in an analogical 
way as described in detail in the recent papers JX0|, [□]]. They are much smaller than in produc- 
tion of light vector mesons. For illustration they are depicted in Figs. |5] and [7| for the lead target 
as a difference between solid and dashed lines at various values of Q 2 . In the photoproduction 
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Figure 7: Nuclear transparency for coherent electroproduction — > 
J/\I/ A as function of energy at Q 2 = 0, 5, 20 and 100 GeV 2 for beryllium, 
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calculations with and without gluon shadowing, respectively. 
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limit Q 2 = the onset of gluon shadowing becomes important at rather high photon energy 
v > 1000 GeV for incoherent and v > 500 GeV for coherent production. This corresponds to 
the claim made in |15| that the onset of gluon shadowing requires smaller XBj than the onset of 
quark shadowing. The reason is that the fluctuations containing gluons are in general heavier 
than the qq and have a shorter CL. 

Although gluon shadowing is included in all calculations, in the kinematic range discussed 
in the present paper it is small enough and consequently does not affect the main achievements 
and conclusions important for investigation of CT effects in coherent and incoherent production 
of charmonia off nuclei. 



6 Summary and conclusions 

In the present paper we focused the main emphasis to production of charmonia due to advan- 
tages as compared with production of light vector mesons [K| . Electroproduction of charmonia 
off nuclei is very convenient for study of an interplay between coherence (shadowing) and for- 
mation (color transparency) effects. A strong inequality l c < If in all kinematic region of v and 
Q 2 leads to a different scenario of mixing of CT and CL effects as compared with light vector 
mesons where l c > // at Q 2 < 1 4- 2 GeV 2 . Consequently, at small and moderate energies a 
problem of CT-CL separation is not so arisen. Besides, due to quite a large mass m c = 1.5 GeV 
the relativistic corrections and nonperturbative effects are sufficiently smaller. Investigation of 
still heavier vector mesons (bottonium, toponium) encounters the problem of very low yields 
as well as very small CT and CL effects (due to very large mass of qq fluctuations and of 
vector meson) to be measured experimentally. Therefore production of charmonia represent 
some compromise because above mentioned theoretical uncertainties typical for light vector 
mesons and very small production rates typical for still heavier vector mesons are eliminated 
to a certain extent keeping sufficiently large CT and CL effects to be detected experimentally. 
This fact supports an enhanced interest to study electroproduction of charmonia off nuclei sep- 
arately. We used from |Kj a rigorous quantum-mechanical approach based on the light-cone 
QCD Green function formalism which naturally incorporates the interference effects of CT and 
CL. Our main results and observations are the following. 



Within suggested approach taken from [10] interpolating between the previously known 



low and high energy limits we studied for the first time CT effects in incoherent and 
coherent electroproduction of charmonia off nuclei. 

As the first test we compare the model predictions with available data from the NMC 
experiment concerning energy dependence of the ratio Tr^/Tr™ of nuclear transparen- 
cies for incoherent production of J/\l/ at Q 2 = 0. We found a good agreement with the 
data what confirms a dominance of CT effects at small and medium and CL effects at 
medium large and large energies. 

The onset of coherence effects (shadowing) can mimic the expected signal of CT in inco- 
herent electroproduction of charmonia at medium large and large energies. In order to 
single out the formation effect data must be taken at such energy and Q 2 which keeps 
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l c = const. Then observation of a rise with Q 2 of nuclear transparency for fixed l c would 
give a signal of color transparency. Predictions of Tr™ c (Q 2 ) as a function of Q 2 at different 
fixed l c shows rather large CT effects in incoherent production of charmonia. Although 
the variation with Q 2 of nuclear transparency at fixed l c is predicted to be less than 
for production of light vector mesons |T(J it is still rather significant to be investigated 



experimentally even in the range Q 2 < 20 GeV . CT effects (the rise with Q 2 of nu- 
clear transparency) are more pronounced at low than at high energies and can be easily 
identified by the planned future experiments. 

The effects of CT in coherent production of charmonia are found to be less pronounced 
similarly as in production of light vector mesons |TI| . A wider range Q 2 < 100 GeV 2 and 



heavy nuclei gives higher chances for experimental investigation of CT. However, on the 
other hand, it faces the problem of low yields at high Q 2 . 

The effects of gluon shadowing were shown to be important only at much higher energies 
than in production of light vector mesons due to large mass of cc fluctuation. Nuclear 
suppression of gluons was calculated within the same LC approach and included in pre- 
dictions. It was manifested that these corrections are quite a small at medium energies 
which are dominant in the process of searching for CT effects. 

Finally one can compare the predictions of charmonium incoherent and coherent produc- 
tion off lead target (see Figs. |], |^, |] and |7|) obtained within rigorous quantum-mechanical 
approach based on the light-cone QCD Green function formalism (incorporating naturally 
CT and CL effects) with results of the paper [|ll[] evaluated in the approximation of long 
coherence length l c ^> Ra (without CT effects) with realistic light-cone wave functions 
of charmonia and with corrections for finite values of l c . We find a nice quantitative 
agreement at moderate and high energies and at low and medium values of Q 2 . This fact 



confirms justification to use that high-energy approximation |TT[ for charmonium elec- 
troproduction off nuclei in the kinematic region where CL effects dominate. Moreover, 
it gives a basis to perform in the future the fully realistic calculations within LC dipole 
approach based on the light-cone Green function formalism employing a realistic dipole 
cross section and realistic LC wave functions of charmonia from [T3| . 



Concluding, the predicted rather large effects of CT in incoherent electroproduction of char- 
monia off nuclei open further possibilities for the search for CT with medium energy electrons 
and can be tested in future experiments. 
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